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Abstract 

We study particles creation from the vacuum by external electric fields, in 

particular, by fields, which are acting for a finite time, in the frame of QED in 

arbitrary space-time dimensions. In all the cases special sets of exact solutions 

of Dirac equation (IN- and OUT- solutions) are constructed. Using them, 

characteristics of the effect are calculated. The time and dimensional analysis 

of the vacuum instability is presented. It is shown that the distributions 

of particles created by quasiconstant electric fields can be written in a form 

which has a thermal character and seams to be universal, i.e. is valid for 

any theory with quasiconstant external fields. Its application, for example, to 

the particles creation in external constant gravitational field reproduces the 

Hawking temperature exactly. 
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I. INTRODUCTION 



The effect of particles creation from vacuum by an external field (vacuum instability 
in an external field) ranks among the most intriguing nonlinear phenomena in quantum 
theory. Its consideration is theoretically important, since it requires one to go beyond the 
scope of the perturbation theory, and its experimental observation would verify the validity 
of the theory in the superstrong field domain. The study of the effect began, in fact, 
first in 3 + 1-dimensional QED in connection with the so-called Klein |I| paradox, which 
revealed the possibility of electron penetration through an arbitrary high barrier formed 
by an external field. Then in 1951 Schwinger found the vacuum-to-vacuum transition 
probability in a constant electric field. It became clear that the effect can actually be 
observed as soon as the external field strength approaches the characteristic value (critical 
field) E c = m 2 c 3 /\e\h ~ 1,3 • 10 16 V/cm. Although a real possibility of creating such fields 
under laboratory conditions does not exist now, these fields can play a role in astrophysics, 
where the characteristic values of electromagnetic fields near and gravitational fields near 
black holes are enormous. One can also mention that Coulomb fields of superheavy nuclei 
can create electron-positron pairs. General considerations, concrete calculations and detailed 
bibliography regarding the vacuum instability in QED can be found in 0-0. Particles 
creation by external gravitational fields and non-Abelian gauge fields can also 

be considered in analogy with electrodynamics. There are also various problems in modern 
quantum theory which are closely related to the vacuum instability in question, for example, 
phase transitions in field theories, the problem of boundary conditions or topology influence 
on the vacuum, the problem of consistent vacuum construction in QCD, string theories, 
multiple particles creation, and so on [|5|,f7t|lOHi5H . 

In spite of the fact that the particles creation effect in external fields was calculated in 
numerous papers, there are still some problems which are interesting to study and discuss. 
In the present paper we are going to focus our attention on the time scenario of the process 
and to consider it in arbitrary dimensions of space-time to be able analyze its dependence 
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on the dimension. To fulfill the first part of the program we consider special external fields 
which act effectively during a finite time and then compare results with ones in a constant 
field. In fact, such a consideration plays also the role of a regularization and helps to solve 
divergence problems which appear in constant external fields. The dimensional analysis may 
be interesting in relation with the study of multidimensional versions of field theories and 
gravity. Lower dimensions, e.g. 2 + 1 dimensions can be of a particular interest. Field 
theoretical models in such dimensions |16|] attracted in the last years a great attention due 



to various reasons: e.g. nontrivial topological properties, and especially the possibility of the 
existence of particles with fractional spins and exotic statistics (anyons), having probably 
applications to fractional Hall effect, high-T c superconductivity and so on [ |17 . 



For calculations we are using the general approach, which was elaborated in the frame- 
work of the field theory for such kind of problems |T8|-f20|,|6| . According to this formulation 
all the information about the processes of particles scattering and creation by an external 
field (in zeroth order with respect to the radiative corrections) can be extracted from special 
complete sets of exact solutions of the relativistic wave equations in the external field (IN- 
and OUT- solutions). A complete collection of exact solutions of such equations in 3 + 1 QED 
is presented in the book ||21|| , in particular, IN- and OUT-solutions and related bibliography 
can be found in ||. That is why in the beginning we analyze and classify exact solutions of 
the Dirac equation in uniform external electric fields in arbitrary space-time dimensions. IN- 
and OUT-solutions are presented explicitly for T-constant, adiabatic, and constant electric 
fields. Probabilities of particles scattering, pairs creation, vacuum-to-vacuum probability 
and mean numbers of particles created are calculated in arbitrary dimensions for three 
types of electric fields mentioned above. The full consideration in the case of the T-constant 
field, which is most important for the time analysis, has no 3 + 1-dimensional analog and is 
presented explicitly for the first time. In spite of some of the formulas in two other cases 
have 3 + 1-dimensional analog, their d-dimensional generalization appears to be non-trivial. 
Moreover, some of these formulas were not presented even in 3 + 1-dimensional case, for 
example, the total mean numbers of particles created and vacuum-to-vacuum probability in 



the adiabatic field. 

Using general expressions obtained for electric fields, which act for a finite time, we study 
the particle creation effect at small and big times. Thus, in particular, we can estimate a 
stabilization time of the process, the time of a pair formation, and give a quasiclassical 
interpretation of a pair creation. Besides, this analysis allows one to select and estimate 
consistently time divergences, which appear in the constant fields. Comparing results ob- 
tained in different time configurations of the electric field, we estimate the role of switching 
on and off effects in the vacuum instability. 

We analyze how the effect of the vacuum instability depends on the space dimensions, 
on the possible boundary conditions, and on a non trivial topology. 

We consider the possibility to add an uniform magnetic field to the electric one and 
calculate the effect. It turns out that one can formulate universal rules to generalize all the 
formulas obtained in the pure electric field to the case when the magnetic field is included 
as well. Its influence on the vacuum instability is studied. 

Finally, it is shown, taking into account the vacuum level shift, that the distributions of 
particles created by the quasiconstant electric fields can be written in a form, which has a 
thermal character and seams to be universal, i.e. is valid for any theory with a quasiconstant 
external fields. Its application, for example, to the particles creation in external constant 
gravitational field, reproduces the Hawking temperature exactly. 

II. GENERAL CONSIDERATION IN AN UNIFORM ELECTRIC FIELD 

The d- dimensional Dirac equation in an external electromagnetic field with potentials 
A^{x) has the form (further h = c = 1) 




(1) 



where ip(x) is a 2'21-component column, 7^ are 7-matrices in d dimensions [22] , 



+ 



V 



= diag(l, -1, -1, . . .), d = D + l 
■> .. ' 
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and x = = (x°, x), x = (x % ), ji = 0, 1, . . . , D, i = 1, . . . , D. 
As usual, it is convenient to present ip(x) in the form 

if>(x) = (P^ + m)^x) . (2) 

Then the functions have to obey the squared Dirac equation in d dimensions, 

(> _ m 2 _ 6 -a^F^ <j>(x) = , F IU , = d,A v - d v A, , = i[ 7 ", 7 1 • (3) 

Let us consider the field F^ u with only one nonzero invariant / = ^F^F^", which sup- 
poses to be negative, I < 0. In this case there exists a reference frame where only the 
components F^ of the field differ from zero. That corresponds to a pure electric field, which 
is a particular case of external fields, violating the vacuum stability (creating particles). Let 
this electric field be uniform. It can be nonstationary, but with a constant direction in the 
space. Then one can always direct it along the Thus, 

F 0l = (0,...,0,E(x )), F lk = 0. (4) 

For such a field we will use the following potentials: A — A\ — . . . — Ad-i — 0, Ad = 
Ad(x°). The constant uniform electric field is of a special interest, because QED with such 
external field (as with any free external field) can be considered as exact QED (without 
external fields) with some special initial states of the electromagnetic field p3| , p4] ,|6[| , which 
provide the corresponding nonzero mean values of the electromagnetic field. Sometimes 
an alternating electric field can also be treated as a slight nonuniform free field, which is 



stipulated by some specific external conditions: existence of a waveguide [25], interference 



of two coherent waves |26[] and so on. However, the study of the constant field shows 



the existence of divergences related to the infinite action time of the field. More correct 
consideration demands a regularization in time. For instance, one can consider a field, 
which acts only a finite time T, being constant within this interval. Such an approach 
allows also to avoid problems with the definition of IN- and OUT-states in non-switching 
external fields at x° — > ±oo. Another possibility is to consider an alternating field, which 



switches on and off adiabatically at x° — > ±00, and is quasiconstant at finite times. In 
the next section we are going to consider all the possibilities mentioned to study the time 
scenario of the particles creation. 

Solutions of the equation (|3]) in the field (^) can be written in the form 

4> P ,sA x ) = 4> P A X °) ex P{^P x }^s,{r}, r = (ri, . . . s = ±1, rj = ±l, (5) 

where f s ,{r} are some constant orthonormal spinors, v s r r \V s ,{r'} — S r>r >. The eq.(||) allows one 
to subject these spinors to some supplementary conditions, 

S±%i,{r} = 0, 3 ± = i(l ± 7 V), rank S ± = J (d) = 2®~\ (6) 

1 X 

R±v s ,{^i,f} = 0, i? ± = -(l±- — -), rankR± = -J (d ), if d > 3, 

2 |Pj_| z 

f = (r 2 , . . . , r^.J, pi = p a , a = 1, . . . , D - 1, = . 

If d < 3 the quantum numbers r do not appear and for d = 2 the perpendicular components 
of the momenta are absent. 

Taking into account the conditions @, one can write an equation for the functions 

<J> P A X °)> 

<M*°) = . (7) 

A formal transition to the spinless case, which corresponds to the use of the Klein-Gordon 
equation instead of the Dirac one, can be done by putting s = in (0) and f s ,{r} = 1 in (U). 

The eq.(^) has two independent solutions at fixed p and s. Thus, an additional quantum 
number ( appears, ( = ± . Combining two independent solutions, which correspond to 
different ( , one can construct two complete sets of solutions ^0 P)S (x°) and ^(j)p tS (x ), obeying 
the following conditions 

i ~d% ^p- 8 ^ ) = c£p e^p.^ )' si § n ( £ p = C> x ° ~* ~°° ' 

i^o C< M X °) = C ^p C 0p, s (x°), sign^ p = C, x°^+oo. (8) 

They provide in turn the following behavior 
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cIxq 



-JI2 + (Pd~ eA D (x )) 2 + pi + m 2 + iseE(x c 



H o .p.{x ) cVw,r(aO = C^p cV'p.s.r^), , sign c £ p = (, x° -> -oo , 

H . p .(x°) c ip p , s , r ( x ) = ( £ p c if> PiS>r (x), sign c £ p = (, x° -> +00 , (9) 

of the corresponding Dirac equation solutions, ^if> P)Str (x) — {jP + m ) t < f>p,a,r{ x ) and 
typ,s,r(x) = (jP + m ) typ,s,r(x)- in the eq. © H . p , = 7°(m + 7P) is one-particle Dirac 
Hamiltonian, and £ are quasi-energies. The solutions ±ip Pj a, r ( x ) describe particle (+) and 
antiparticle (-) in the initial time-instant whereas ± ip PyS ,r{ x ) describe particle (+) and an- 
tiparticle (-) in the final time-instant [19, Gj. 



One can see that the solutions with different s and fixed (, p, r are dependent. For 
example, 

(1pp,s,r(x) = m + C^p- S ,r(z), (10) 

( a p,-s 

where b Pyr = ri|pj_| if d > 3, b Pir — p\ if d — 3, b Pir = if d = 2, and ^a PjS are some 
coefficients. To see how fllOD appears one can use (||), (^), @ and the following consequence 
of two latter, 



i-j^ + s(>d - eA D (x )) 



(<PpA x °) = < a P,s t<t>p,-s{x°) , c a p _i c a Pi+ i = m 2 + p 2 ± . 



Similar relation holds for ^if} p>sr (x). The eq . C|X0|) means, in fact, that the spin projections of 
a particle (+) and an antiparticle (-) can take on only Ju\ values. Taking that into account, 
one can only use the following independent solutions, 

±%,r{ x ) = ( / y p + m )±0 P ,±i,r(z) , ± ^p,r( X ) = (l P + m ) ± P ,Ti,r(a ; ) • (11) 

Further, we are going to calculate different matrix elements between the solutions (|TT]) by 
means of the conventional time independent Dirac scalar product (ifj, if}') = J ^{x)j if}'(x)dx. 
In the case under consideration, due to the above mentioned properties ©,(0) of the spinors 
C</>p,s,r(aO and ^<p PjStr (x), the scalar product can be reduced to a form which is convenient for 
calculation and, in particular, does not contain 7 matrices at all, 

(+i>p,r,+if>p>,r') = i(2ir) D 5 r y5(p - p) ;0 P>+ i(x°) d (ido+PD - eA D (x )) ~0 Pi+ i(x°), 
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(tW,lVy,rO = t(2n) D 5 r , r/ 5(p - p') U P ^(x ) X (id - p D + eA D (x )) ^0p,-l(x°), 

(T^p.r, t^y) = l(27T) D 6 r y6(p - P ')(m - ib p>r ) l<f> P>+1 {x°) X ±<j> P) -l(x°) , (12) 



where c? = "9o ~ <9o • (The right side of fll2"|) reproduces the corresponding Klein-Gordon 
scalar product if one puts formally C[^o ± (pd — gAd(x ))] — m — ib — 1.) 

One can see from (|) and (|P2"D that the solutions (|TTJ) can be normalized to obey the 
orthonormality relations, 

(cVV> C'^p'.r') = S C ^5 r y5{p - p'), ( C ^ P ,r C V P ',r') = 8 r ,r> 8 (P ~ p')- ( 13 ) 

Moreover, each set of solutions ^ip Pjr (x) and ^^) pr (x) forms a complete system, thus, we are 
dealing with the so-called IN- and OUT-sets of solutions correspondingly |T9|j20| J^]. 



Using ([TJD, one can find decomposition coefficients G [cl^'j °f the OUT-solutions in the 
IN-solutions, 

ty(x) = +^(x)G( + | c ) +_^(x)G(_| c ) . (14) 
The matrices G fcF') ODe y the following relations, 

G(c\ + )G( C \+V + K G( < \-)G( < \-V = (0)^, 



G ( + | + ) G ( + | + ) t + ( + |-) G (_|") t = , (15) 

where I is the unit matrix and k = ±1 for fermions and bosons respectively. Relations ( |T5D 
can be derived from the conditions (13). Due to eq. flT^) we can easily see that the matrices 



G \ c\ < ") are diagonal, 



All the information about the processes of particles creation, annihilation, and scattering in 
an external field (without radiative corrections) one can extract from the matrices G 
because they define a canonical transformation between IN and OUT creation and annihi 
lation operators in the generalized Furry representation |19|,20,6||, 



a! (out) = a ] {in)G ( + |+) + b(in)G 
b(out) = a\in)G (+|~) + b{in)G (_|~) . 



(17) 



Here a^in), b^iin), a n (in), b n (in) are creation and annihilation operators of IN-particles and 
antiparticles respectively and a^out) fi^out) , a n (out), b n (out) are ones of OUT-particles 
and antiparticles, n presents momentum p and spin projections r. For example, the mean 
numbers of particles created (which are also equal to the numbers of pairs created) by the 
external field from the IN- vacuum |0, in > with a given momentum p and spin projections 



r is 



N Pi r =< 0,in\a pr (out)a Pir (out)\0,in >= g (-| + ) 



(18) 



Here the standard volume regularization was used, so that 6(p — p') — > S pp '. The probabil- 
ities of a particle scattering and of a pair creation have the following forms respectively 



-P(+I+W,p',r' — I < 0, out\a p ^.(out)a} , r ,(in)\0, in > | 2 = S r}r '5 P}P 



1 



1 - kK 



-P, 



(19) 



p.r 



N, 



P{— + |0)p,r, P ',r' = | < 0,out\bp tr (out)a p > >r '(out)\0,in > I = 5 rtr >8p tP >- '—— P v , (20) 

where 10, out > is the OUT vacuum and 



Pv 



< 



0, out\0, in > | 2 = exp |«^ln (1 — /tiVp )r ) j 



(21) 



is the probability for a vacuum to remain a vacuum. The probabilities for an antiparticle 
scattering and a pair annihilation are described by the same expressions P(+|+) and P( — h 
|0) respectively. 

Thus, to be able to calculate the quantities (|T8|)-(|2T|), in the case under consideration, 
one has to find solutions of the ordinary differential equation (|7|), which is in fact Schrodinger 
equation for a linear oscillator with time-dependent frequency. However, one can make some 
general conclusions, which do not depend on the concrete time dependence of the electric 
field in eq.(|7|). First of all, the matrices G N^j are diagonal in all the quantum numbers 
introduced. Second, the quantum numbers r do not enter in the eq.(0) and, due to the 



structure of the scalar product (0), the matrices can depend on r via a phase only. That is 
why all the probabilities and the mean numbers do not depend on r, so that in the fermionic 
case the total (summed over all r) probabilities and the mean numbers are Jw) times greater 
than the correspondent differential quantities. For example, the total number of particles 
created with a given momentum p is 



^2 Np,r = J(d)Np !r . 



(22) 



Finally, it is clear that due to the structure of the eq. (|7|) and the scalar product fli~2]) the 
dimensionality d enter in the differential probabilities and mean values via the combination 
Pi only. 

III. T-CONSTANT, ADIABATIC, AND CONSTANT ELECTRIC FIELDS 

A. T-constant field 



To analyze the time dependence of the particles creation effects let us consider the field 
(|j) with E(x°) having the form 

o, x° e I 

E(x°) = \ E, x°ell (23) 
0, a; e in , 

where the time intervals are: I = (—oo,ti), II = [ti , t 2 ] , HI = (h, +oo), t 2 — ti = T, t 2 = 
—ti, and eE > is chosen. Thus, in fact, we consider a constant electric field E, which is 
acting a finite time T. Further we will call it T-constant field. The corresponding potential 
Ad{x°) can be chosen in the form 



A D (x°) = 



Et u x° el 
Ex , x° e II 
Et 2 , x° e III . 



(24) 
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In each interval I, II, III the equation (0) has two independent solutions, which are 
correspondingly in J: exp{— ip (ti)x } and expj+zpo^^x }, in 77: D v _i±e [(1 — and 
D_ v _i^b_ [(1 + and in J/J: exp{— ipo^)^ } and exp{+2p (^2)^ }, where D v (z) are We- 
ber parabolic cylinder functions (WPC-f unctions) [p^j , and 

1/ = y, A = — — — , f(ac ) = — — , po{x ) = \/m 2 + pl + (p D - eA D (x )y . 

Using them and conditions one can construct IN- and OUT-solutions -ip Ptr (x) and 
+ ?/'p ir (a;) (see Sect. II). The corresponding expressions for _</) p> _i(x°) and + Pi _i(x°) are of 
the form 



exp{+z*po(^i)(z° - h)}, I 
Ci { ai D u [(1 - z)£] + 02£)_ v _i [(1 + z)£] , // 
</( + |_) exp{-^ (t2)(x° - t 2 )} + g(-\-) exp{+tp (t 2 )(x° - t 2 )}, /// ; 



(25) 



(26) 



>p,-i(^°) 

</( + |+) exp{-zp (ti)(x° - tx)} + <?(_|+) exp{+ip (ti){x° - h)}, I 
C 2 < a' x D v [(1 - i)£] + a' 2 D^ [(1 + z)£] , II 
exp{-ip (t 2 )(x° - t 2 )} , III , 

where the normalization constants are 

d = (2n)- D ^(2p (t t )q t )- 1 / 2 , Qi = p (t t ) - (-l)\p D - eA D (U)), i = 1, 2. 

To provide the continuity of the solutions in the time instants t\ and t 2 one has to impose 
the following conditions 

+0 p ,_i(t 4 + 0) = j# p ,_i(ti - 0), —t^iu + 0) = —t^ x {u - 0) , 

which allow one to define step by step all the coefficients a«, a[, and g (±| + ), g ( ± |-) • The 
first ones are 



a, = -(-V 



MV2eE 



Mt2)fj-\t 2 ) 
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where 



n n 1 77 

M = D v {z) —D_ p _i(iz) - D_ v _ x {iz)—D v {z) = exp{-(z, + 1)-} 



is the Wronskian determinant 127 , and 



/IV) = |i± 



id 



D-„-x [(i + m , / 2 (±) (^°) = ( i ± -% ) d u [(i - m 



They can be used to define the latter coefficients. From those we need to know explicitly 
only g (_| + ) and g ( + |-), which are 



9 - 



exp{(z/ + 1) — } 



i f Po(,ti)qip {t 2 y 



8eEq 2 



/ a H (*2)/i H (*i)-/r(t2)/r(*i) 



(-) uHi 



exp{(, + l)^}f Po(tl)92Po(t2) ' 
PU ; 2 1 \ 8eE qi 



/ 2 (+) (t 2 )/i (+) (ti)-/r(t 2 )/r j (ti) • (27) 



One can see that the coefficients (|27| ) obey the properties 



Pd->~Pd 



(28) 



The first one can be verified directly, whereas the second one is easy to derive comparing 
representations of the scalar product (|12"D in the time instants t\ and t 2 . Thus, one can 
conclude that \g (-| + )| is an even function of the momentum pr>. 

To calculate the probabilities and the mean numbers according to the formulas ( [T8| - |2lD 
we need really to know only the coefficients Comparing (P^) and (p5|), (|26D, we 

conclude that the relations ( p!6|) hold and we have the expression for the mean numbers 
of pairs created in the form flTSP, in which g (_| + ) is defined by (P7|). In fact, the mean 
numbers N p>r define all the probabilities via the formulas (|l9l - ^Tl) . As it was shown above 
this function is even in all the momenta p, including p^ and does not depend on the spin 
quantum number r. Using the recipe presented in Sect. II, it is easy to get an explicit form 
iVp for the the bosonic case from the fermionic one (|18D and (^7|). 

Now we are going to analyze the dependence of all the characteristics on the time T and 
on the momenta. One has to remark that the dependence of the longitudinal momentum 
Pr> is of a special interest. This dependence is essentially correlated with the T-dependence. 
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One can see, e.g. from the four-dimensional case [I8|,6|, that in the constant field (T = oo) 
all the characteristics do not depend on the momentum pp. This is a source of some kind 
divergences if one is interested in the total characteristics, which require integration over po- 
Due to the reasons mentioned above it is enough to analyze only the quantity N Ptr . As to 
the momentum p^, one can restrict itself only by pr> positive or p£> negative. Remark that 
the momenta po enter in all the formulas via two dimensionless parameters £1 and £2 only, 



which in turn appear in the WPC-functions. Thus, in fact, one needs to analyze the depen- 



dence of the latter on £. It is convenient to consider the region —yeE^ < £1 < +°° > £2 > 



eE~, which corresponds to po negative, < — po < +00, in particular, there always 
In the region £1 > K, where if is a given number K » 1 + A (in terms of the momentum 



this region corresponds to \p^\ > eEjr + KyeE) one can use the asymptotic expansion of 



WPC- functions |7|, 

A,(z) = ^exp{-^/4}(£ 



N 

c 

n=0 



n\(-\z 2 ) n 



to conclude that at any T the behavior of the mean numbers ([18] 

3 



argz| < -7T, 



is 



(29) 



N p , r = O 



(30) 



For small T, T << 



,,, and \pd\ « eE^, one can get 



eET 2 



1 + 0(v / e£T) 



(31) 



p ' r 4A + eET 2 

At T << and T << the form flJH]) reduces to iV Pjr = eET 2 /(AX) and coincides 

with one which can be derived in the frame of perturbation theory with respect to the 
external field. 

The most important region for the time divergences is one of big T, namely, let us 
consider T >> ^7=^(1 + A). In this case £ 2 is always big and positive £ 2 >> 1 + A, so that the 
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asymptotic expansion ( p9]) can be used for any given momentum pr>. As to the parameter 



£1, the whole interval —yeE^ < £1 < +00 can be divided in three regions: 



a) - VeE- < 6 < -K\ b) - K < & < K; c) & > K . 
The mean numbers N p>r were estimated in the region c) before, see Q50p. In the region 



a) one can use some relations between the WPC-functions |27j, for example, D u (z) = 
exp{iTcis}D u (— z) — f(Z^y exp iz), and the asymptotic expansion (|29f) . Then one 

finds 



N = e 



1 + 



1 + A 
6 



+ 



l + A 

6 



' < £1 < (32) 



The latter expression allows one to consider the limit T — > 00 at any given p. In this 
limit the mean numbers take a simple form 



N = e 



-7TA 



(33) 



Thus, when the electric field is acting for a long enough time, the mean numbers of particles 
created in a given quantum state are stabilized and coincide with expressions which were 
obtained in the constant electric field in 3+1 QED ]Tj|. (The stabilization (^) was first 
remarked in f28| for particles created with zero momenta, using a finite action electric field 
in 3+1 QED.) One can also estimate a characteristic time of such a stabilization. To this end 
one can see that is a small parameter in the decomposition (|32|) in case \p D \ « eE~. 
If T » T , T = ^± then the mean numbers are stabilized and T is the above mentioned 
characteristic time. 

The intermediate region b) does not allows one to use an asymptotic expansion of WPC- 
functions to analyze the £1 dependence. However, one can make some conclusions about its 
contribution in integrals over the momenta. For example, due to the Fermi statistics N v ^ r is 
always smaller then unity, that is why the integral over the momentum pr> in the region b) 
is less then 2\feEK and is not essential in comparison with the same integral in the region 
a) at T — » 00. 
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Using these considerations, one can now estimate the sum over the longitudinal momen- 
tum pd of N p>r , which is the mean number of particles created with all possible momenta 
Pd- To do this we go over to the integral, J2 PD "~ > jr I dpo, where L is the length in the 
direction x D . As was shown above, at T >> ^7=g(l + A), N p>r is quasiconstant in the area 
a), the asymptotics in the area c) has the form ( |30D and the contribution to the integral 
form the area b) is less then ^/eEKL/ir. Then one can conclude 



L r +oc veEL r / — x n 
N p± , r = — / N p , r dp D = \ \feETe + 0{K) . (34) 

Thus, when T >> -j=§K » -^(1 + A) we can effectively replace the integral over po by 
eET and write 

N p±tr = A ion9 e-" A , A lon9 = ^-eELT . (35) 

Z7T 

The factor Ai ong can be can interpreted as the total number of states with the longitudinal 
momenta pr, of particles created. 

It turns out that the expressions for N p ^. and Ap ± r at big T for scalar particles coincide 
with the ones for spinor particles. 

To get the total number A^ of particles created one can sum over the spin projections, 
using eq. fl22|) , and then over the transversal momenta. The latter sum can be easily trans- 
formed into an integral, 

^ = E^p=(^t/^p, (36) 
where V(d-i) is (d-l)-dimensional spatial volume. Thus, on gets 

where E c = m 2 / e is the critical field strength. As one can see the velocity of particles creation 
is constant at big T. The same result was obtained in 3+1 dimensions in the framework of 



the QED at finite temperature, using the functional Schrodinger picture p9[ . 



15 



The vacuum-to- vacuum transition probability fl2T|) can be calculated, using both kinds of 
regularizations, with respect to the volume and to the time. Thus, we get the d-dimensional 
analog of the well-known Schwinger formula [f2| , 

-l^f r k, 

n=0 (71 + 1)2 

Returning to the Schwinger result, one ought to say that it was, in fact, obtained from the 
constant field consideration by means of a regularization. Thus, the space-time volume VT 



00 (— l)( 1_re )f < E c \ 
P v = exp {-/iiV} , ^ = E —t- ex P \ - nn -fF f ' ( 38 ) 



appeared in his formula. One ought to mention the paper |15[ where the Schwinger result 



in the constant electric field was recovered, using the functional Schrodinger equation. 

B. Adiabatic field 

Let us consider an alternating uniform electric field (|]), where the function E(x°) has 
the following form 

E{x°) = £cosh- 2 ^ . (39) 

Such a field switches on and off adiabatically at x° — > ±00, and is quasiconstant at finite 
times. We will call it adiabatic field. The corresponding nonzero potential is 

x° 

A D (x°) = a£tanh— . (40) 
a 

In this case solutions of eq.(0) can be written in terms of hypergeometric functions 
F(a, b; c; y) [27], for example, 

c p , +1 (x°) = c Ce---*° (l + e^J 2 ~ + c u(x°), (41) 
+u{x°) = F(a, b; c; y), -u{x°) = y l ~ c F{a - c + 1, b - c + 1; 2 - c; y), 
a = — (2eEa + uj + — uo-), b = 1 + — (— 2eEa + uj + — cj_), 



1 x° 1 

c = 1 — zau;_, y = -(l + tanh — ), u>± = dm 2 + + (pn T eEa) 2 , 

2 a v 

where are some normalization constants. Considering the asymptotic of the functions ( f4"T| ) 
at x° — > —00, (in this case F(a, b; c; y) — 1 f27|), one can verify that the relations (f|) hold and 
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(£ p = C^-- Moreover, the solutions at x° — > — oo describe free particles. By analogy one 
can construct solutions ^<p P)+ i(x°), which describe free particles with energies ^£ p = (u + 
at x° — > +00. Their asymptotic form at x° — > +00 is ^0p ;+ i(x o ) = ^Cexp(— iQuj + x Q ). 
To calculate the coefficients G (+\~) from (0) it is enough to know the corresponding 
asymptotic + P;+ i(x°) and ~ <f) Pt+ i(x°) , let say at x° — > +00, and normalization constants 
"C = (27r)- D / 2 (2cu + (cu+-p D + eE«))- 1/2 , = (2tt)- d / 2 (2cu_(cu_ + p D + eEa))~ 1/2 . 

Thus, we get the mean numbers of fermions created, 



sinh 



(2eEa + 



uj 4 



sinh 



7TM 

2 



(2e£a + 



p.r 



(42) 



sinh(7rao; + ) sinh(7rau;_) 
In 3+1 QED the corresponding formula was found first in |3~0|| . For scalar particles it has a 
different form 



cosh 



71 J (eEa 



2\2 _ I 

4 



+ sinh 



p.r 



(43) 



sinh(7rau; + ) sinh(7rau;_ 
Let us consider the a-dependence of these expressions. For small a, a « ^i/m 2 + p 2 , 
when the potential changes sharply, we get for fermions 

(neEa 2 ) 2 [ I 



N, 



m 2 +p 2 



p.? 



sinh (71a J m? + p 2 ) 



(44) 



and for bosons 



{neEa 



2\2 



{eEa 2 ) 2 + 



m 2 +p 2 



p,r 



sinh (naJm 2 + p 2 ) 



(45) 



Small a in the case under consideration corresponds in a sense to small T of the T-constant 
field. Thus, we have to compare the expressions ( pi]) and ([44]). One can see that they are 
quite different, so that the effects of switching on and off are essential at small times. 

Further let us consider big a only, a » ^r=g(l + VA). Then the mean numbers for 
fermions and bosons have the same form, 



N Ptr = exp {— na(uj + + uj_ — 2eEa)} 



Let us take small longitudinal momenta \pn\ << eEa, then 



(46) 
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N pr = exp < — 7rA 



1 



KeEaJ 



(47) 



Considering the limit a — > oo, one gets the formula fl33|). That means that the effects of 
switching on and off are not essential at big times and small longitudinal momenta. For big 
longitudinal momenta \po\ » eEa, the mean numbers of particles created are exponentially 
small, 

iV Pir = exp {-2ira(\p D \ - eEa)} . (48) 

Let us find the total numbers of particles created with all the longitudinal momenta at any 
fixed p_i_, r. Passing from the summation over pp to the corresponding integration, we get 

Comparison with the formula (j35|) shows that the adiabatic field at big times (big a, a » 
7=(1 + vT) and fixed p_i_, r is equivalent to the T-constant field at T >> -r*g(l + A) with 



the identification a = y/\T. To do summation over all transversal momenta, it is convenient 
to use the representation 

1 f°° 
—j= = 21 exp(—7r\s 2 )ds. 
VA Jo 

Then the total number N of the particles created reads 

i±^ V {d _ 1) a5m d (E\i ( E c \ 
N = J (d) TdwTl hFT exp -7T— , (50) 



'« (2^p \Ej 

where 

/•oo t d _ 2 ^ / i 

5=/ dtri(t + l)-— exp(-tvr-^) = v 7 ^^ - 
jo -£/ \ 2 



^_ /l d-2 E c 

: 7r — 

2 ' £ 



is expressed via the confluent hypergeometric function ^7|. The vacuum-to- vacuum transi- 
tion probability P v has the form 



P„ = exp {-/iiV} , // = ^ -P— exp ^ -7wr-^ [ , (51) 

n =o (n + 1) 2 L & 



c -i / — t i d - 2 ^ 
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If E/E c << 1 one can use an asymptotic of \I/-function f27j, \I/ (|, — mr^j = -j=^Jf + 

3/2 N 



o 



Then 5 = yjs e n = n~2 and /x = 1. In this case the adiabatic field is 
equivalent to T-constant field with the identification a = T^J~^-. At strong fields E ~ E c 
all the terms with different e n contribute to the sum in ([)]]) and the expression for P v differs 
essentially from the one for the T-constant field. 

One can also remark that the case of a periodic alternating electric field in 3+1 dimensions 
was also considered in the literature, for example, quasiclassically [pll], and exactly [32 . 



C. Constant field 

Here we consider the case of a constant uniform electric field @). Then E(x°) = E and 
potential = Ex . Solutions of the eq.(^) in such a field can be found in the form 

;0 p , s (x°) = CD V _^.(±(1 - i0 P >°) = CD_ u _^.(±(l + i)^). (52) 

Using an asymptotic expansion of WPC-functions fl29|) , one can get the asymptotics of the 
quasi-energies, 

( S P = (\eEx° - p D \, ? £ p = ((eEx° - p D ) , 

so that IN and OUT-solutions can be constructed from ( |52]) by means of eq. QTT|). The 
same asymptotic expansion ( |2^D allows one to calculate the normalization constants, C = 
(2n)- D / 2 (2eE)- 1 / 2 exp {-ttA/8} for spinor case and C = (2n)~ D / 2 (2eE)- 1 / 4 exp {-ttA/8} 
for scalar one. Straightforward calculations, similar to ones where made in the two previous 
cases, lead to the expression (|33|) for the mean numbers of particles created. It does not 
depend on the dimensionality of the space and coincides with the result which was derived in 
|T8| for 3+1 QED. In that paper the authors used quasiclassical considerations to advocate 
the classification of the solutions (|52|). The constant character of the field does not allow 
one to treat consistently time divergences, so that they got over them "by hand" , using also 
quasiclassical considerations. Now one can see that the consideration of the T-constant field 
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gives a possibility both to ground all the results obtained from the constant field solutions, 
solving consistently the problem of the time divergences and to go beyond the scope of the 
constant field to analyze the time scenario of the process. 



D. Inclusion of a magnetic field 

In the same manner as before one can consider a more general case when a constant 
uniform magnetic field is included, provided the invariant I is negative. In fact, in d > 3 
there are [|] — 1 independent invariant parameters Hj, j = 1,2, ... , H^_ 1 of the magnetic 
field, that corresponds to the possibility to construct [~] invariants of the electromagnetic 
field. In a convenient reference frame, the magnetic part of the field tensor F^ v is presented 
by the components F£, = Ej=i Hjitf+Hi - One can always select solutions of 

the squared Dirac equation (^|) as eigenfunctions for all independent nonzero terms, which 
describe the interaction of intrinsic magnetic moment of a particle with the external magnetic 
field. In this case the matrices G f^'j are diagonal and one can construct them using the 
corresponding expressions in the pure electric field. Namely, for d > 3 one has to make there 
a replacement 

t^" 1 0, d is even 

i=i p 2 d _ 2 , d is odd , 



\eHj | (2nj + 1 - r,-), % = 0, 1, . . . , H i 7^ 
Pj+Pj+i > #i = ° 

In the presence of the magnetic field some momenta pj have to be replaced by the discrete 
quantum numbers rij. The number of these momenta Pj corresponds to the number of 
nonzero parameters Hj. The magnetic field lifts the degeneracy in spin projections in all 
the characteristics of the particles creation effect. 

We present here explicit formulas in presence of the magnetic field for the total charac- 
teristics iV and P v in case of the T-constant field at big T, 
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7V " J W 2(^-1)^/2 ^ eXP l E 



J c 
oo 



P. = exp { -,N } , /i=E «P {—§ J , (54) 



where 



= g [ sinHnnk/E] ^^/E)] - j , d is even , 

/ m 2 E\^ {d ~ 3)/2 { eH 
/3(n) = — - TT ^ - / J rr /rnN [coshfnTrF./S)]— , d is odd . 

The corresponding formulas for 3 + 1- dimensional case were first written in |18,3|, and, in 



fact, can be derived easily from the calculations of Schwinger 0. They follow from (54) at 
d = A. 

There exists a possibility to get also exact results for the particles creation in case when 
a plane wave is added to the combination of electric and magnetic fields. The corresponding 
calculations for (3+l)-dimensional case were made in cited6 for the electric field plus a plane 
wave field, and in |34j] for a general combination of electric, magnetic, and plane wave fields. 
They can be generalized to any dimensions, combining the approaches of |34[ and of the 
present paper. 

IV. DISCUSSION 
A. Time and space-dimensional analysis 

The calculation and analysis presented in Sect. Ill for fields, which are effectively acting 
for a finite time, allows one to study both pairs formation in time and the role of switching 
on and off effects. Besides, due to the fact that these calculations are made in arbitrary 
dimensions of the Minkowski space-time, one gets a possibility to analyze the influence of 
the dimensionality on the vacuum instability. 

Studying the T-constant field, one can see that the stabilization of the mean numbers of 



particles created with given p, r in the form (]33|) for the longitudinal momenta \pd\ << e E 



2 
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comes at T >> To, where To = + A). The characteristic time T can be called 

stabilization time. At the same time N p>r for the big longitudinal momenta \po\ » e-E 1 ^ 
decrease according to the power low (|30|). 

The stabilization of the mean numbers with given p, r in the adiabatic field in the same 
form (|33|) comes for the longitudinal momenta \po\ << eEa at a » cto, cto = ^=(l + vA). 
For big \pr>\ » eEa the mean numbers are exponentially small fl48f) . For big a the adiabatic 
field varies slowly and coincides nearly with the constant one in the time interval < a. 
Then «o is a characteristic time of the stabilization in this field. Thus, the stabilization 
time «o i n the adiabatic field differs from the corresponding time To in the T-constant field. 
Thus, one can believe that the stabilization process depends of the switching on and off 
effects. In the case E/E c < 1, which corresponds to A > 1, one can see the stabilization 
comes quicker for adiabatic field than for the T-constant one (a^ < T ), i.e. the adiabatic 
form of switching on and off affects less the quantum system than the instantaneous one in 
the T-constant field. If E/E c > 1 , there exists a domain of the transversal momenta 
where A < 1. In this case the stabilization times in both cases are the same, olq ~ Tq ~ -7^7?, 
so that for any E the relation a < T holds. 

Thus, one can conclude, that in some cases T-constant and adiabatic electric fields act 
on the vacuum in a similar way. However, the momentum dependence of the mean numbers 
N Pt j. differs essentially at big momenta for both fields. That is related to the switching on 
and off effects. To estimate the role of the effects of switching on and off on the whole it is 
convenient to compare total characteristics. First, let us compare the total mean numbers 
with all the longitudinal momenta, namely, compare the formulas (j35| ) and (f49[). In this case 
the effective action of both kind of fields is the same if to identify a with a/AT. In spite of 
this identification of a and T is different for different A (for different pjj, one can use it in 
a domain pj_ of the transversal momentum, |Ap_|_| << yVn 2 + p^_. As to the total numbers 
(|37|) and (|50|), they coincide if one accepts the identification a = TS^ 1 . However, this 
identification provides only the coincidence of P v for both cases ( |5%| ) and ( |oT| ) if E/E c << 1 



22 



(then 5 = yjf-)- In this case the coefficients /i in (^) and (|5l"l) are the same. One can 
conclude that the effects of switching on and off are not essential for E/E c << 1 and for 
big T >> ^ (-ff) ^ j or for big a respectively. In case of strong fields, E/E c > 1, these 
effects appear to be essential and one has to take into account the back reaction of particles 
created for more realistic external field definition (e.g. see [^j|7|,^4|). 

The stabilization of the mean numbers of particles created with given p,r at T >> T 
can be interpret in the following way: In the T-constant electric field at the finite time 
instant ^ the pairs are created with quantum numbers \pd\ < eE^ (equal for particles and 
antiparticles). This corresponds to the region < 7r D (^) < eET of the observed kinetic 
momenta, 7r°(?) = — (pr> — eAo(~)) = — po + eE^ of a particle in each pair (direction of 
antiparticle kinetic momenta is opposite). At T >> T the effects of switching on and off 
are already not essential. That is why the probabilities of pairs creation do not depend on 
the time instants One can think that at this time instant the particles in 

the pairs are materialized with almost zero longitudinal kinetic momenta at any given p^, 
i.e. with the energies \Jm? + pj_. Then the electric field accelerates them until the end of its 
action. Let us suppose that a particle was created at a time instant t. An expression for the 
kinetic momentum of such a particle at the final time instant (which is equal to its expression 
in the time instant when the field switches off) can be found solving the classical equation of 
motion ^jj- = eE, so that 7r D (|-) = eE — t). Thus, a particle, which was discovered with 
the quantum number po at the final time instant, was created at the time instant t — ^B. 
Then the integration over the longitudinal momenta pr, is equivalent to one over the time 
t, J dpo = eET. This conclusion coincides with one derived in course of the strict quantum 
analysis presented in Subsect. IIIA. According to the same interpretation, for particles 
with relatively nonzero mean numbers, the maximum value of the kinetic momentum 7r D (-|) 
is eET that corresponds to the particles, which were created at the initial time t — — ^, 
whereas its minimal value is and corresponds to the particles, which were created at the 
final time instant t = This conclusion coincides also with one derived from quantum 
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consideration in Subsect. Ill A. 

In the conclusion of the time analysis, one can remark that the time T , which was 
introduced by us as the stabilization time, was interpreted in some papers as the time of 



time «o is different, thus T is not an universal characteristic and depends of the field form. 
In this connection one can propose another characteristic time, which a pair formation in 
a quasiconstant electric field. Indeed, as we have mentioned above, all the results in the 
T-constant and adiabatic fields are comparable if E/E c « 1. In this case the adiabatic 
form of the field is disturbing the quantum system less than the T-constant field. Here 
«o << T and «o ~ Tq = . Since the adiabatic field is quasiconstant for the time 
interval Tq and it is big enough for the stabilization, one can interpret Tq as the time of a 
pair formation. One can extrapolate this interpretation of Tq for any field strength E. A 
quasiclassical consideration confirms this interpretation. Thus, a virtual particle with initial 



zero energy gets from the electric field for the time Tq the energy ym 2 + p^ necessary 
for the materialization. It is easily to see that the time Tq is always either less than the 
stabilization times To, ao or equal to them. Some other consideration related to the time Tq 
one can meet in the next subsection. 

Turning to the dimensional analysis, one can see that the increase of degrees of freedom, 
due to the increase of the dimensionality of the space-time itself and due to the related 
increase of the spinning space dimension J^), affects essentially the total numbers of particles 
created in the unit of the volume and the probability for vacuum to remain a vacuum. Thus, 
the increase of spinning degrees of freedom leads to an increase of N and P v at any ratio 
E/E c . In particular, in d > 3 the numbers of fermions created is greater than the one 
of bosons. An increase of spatial dimensions leads to a decrease of the total numbers of 
particles created in the unit of the volume and the probability for a vacuum to remain a 
vacuum at E/E c < 1 and their increase at E/E c > 1 . 

The presence of walls or of a nontrivial topology affects the spectrum of particles cre- 
ated. If the length Li of the space in the direction of an axis x % is restricted by the walls 



a pair creation p8|j28[] . However, we have seen that in the adiabatic field the stabilization 
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the corresponding momentum \pi\ = 2h*, n = 1,2, . . . is quantized. At Lj ~ — the de- 
pendence of the mean numbers on the boundary conditions is essential. At Li « -j= the 
mean numbers N p%r in the quasistationary fields are very small for any strength E. In this 
connection one can treat Lq = — [1 + (-^) 1//2 ] as a characteristic dimension of the system, for 
which the boundary conditions are essential. At E/E c > 1 it is the Compton wave length. 
It is interesting to remark that the stabilization times T , «o coincide with L Q at E/E c = 1. 

Imposing periodic conditions in the direction of an axis x l (that corresponds, in par- 
ticular, to the torus topology), one gets for the momentum \pi\ = n = 0, 1,2, . . . . 
Then at Lj << ^7= only particles with pi = can be created. If the electric field has the 
same direction, then the total number N and the probability P v do not depend on time T, 
since this dependence arises in course of a summation over the longitudinal momenta. It is 
interesting that the presence of the magnetic field acts as a dimensional reduction. Indeed, 
in the strong magnetic field with some Hj » E the lowest energy level of a boson can not 
be less than \eHj\, whereas for a fermion it can. That means that the strong magnetic field 
acts on bosons as some walls and on fermions as the presence of the torus topology. Thus, 
one can see that if some of the magnetic fields are strong enough, then the corresponding 
spin projection become frozen and total characteristics, like total mean numbers decrease. 
These dimensional effects may be relevant to the matter creation at early universe. 

B. Relation between the vacuum instability in external electromagnetic and 

gravitational fields 

It is interesting to compare particles creation in external electromagnetic fields and in 
external fields of different nature, for example, in external gravitational fields. To this end 
one can use results obtained in the quasiconstant electric fields and in the static gravitational 
fields. The latter problem was considered first by Hawking || who, in particular, calculated 
the mean numbers of particles created by static gravitational field of a black hole with mass 
M in a specific thermal environment, 
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9 U 

exp < 2ir- 



(55) 



9(H) J 

where ui is the energy of a particle created, which supposes to be dependent on a complete set 
of quantum numbers n, g^n) = ^r, where r g is the gravitational radius, so that g^u) is free 
falling acceleration at this radius. This spectrum was interpreted as a Planck distribution 
with the temperature 6(h) — J~jp~ (&b is the Boltzmann constant). As before k — +1 for 



fermions and k = — 1 for bosons. It is also known |35j that an observer, which is moving with 
a constant acceleration gm\ (with respect to its proper time), will register in the Minkowski 
vacuum some particles (Rindler particles). The mean numbers of Rindler bosons have the 
same Planck form (|55|) (with k = —1), where one has to replace g^) by (?(_r), so that 
the correspondent temperature is = ^f^- O ne can ^ n( ^ m any other examples when 
the particles creation in external gravitation fields (and due to a nontrivial topology) can 
be described by means of an effective temperature ]5],[7| ( see also references in the recent 
publications |36[| ). On the other hand the distributions obtained in external electromagnetic 
fields have not the thermal form at a first glance. Nevertheless, there were attempts to 
find close relations between the distributions in both cases, moreover, to derive Hawking 
distribution from the one in external electromagnetic field. 



In |3jJ , the distribution (|33|) at pj_ = was interpreted as the Boltzmann one for particles 
in the ground state with the energy m and the effective temperature 6(e) = The 



same temperature follows from some other consideration |$8| for the same restricted case. 
Unfortunately, such an interpretation does not allow one to include other states with nonzero 
momenta in the consideration. 



In [39], the authors did not introduce an effective temperature directly in the electro- 
dynamical case but tried to find a relation between both distributions, in particular, to 
extract the Hawking temperature from the electrodynamical formulas. We are going to re- 
peat briefly here this consideration , using some new details, which came from the results 
of the present paper. As was established, a particle with given momenta is created in a 



time instant with energy uj = ym 2 + p\, which corresponds to zero longitudinal kinetic 
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momentum at this time instant. Thus, namely this expression plays the role of the total 
energy of the particle at the time moment of creation. Then we can compare equations of 
motion for a classical particle in a constant electric field dn/dx = eE with ones in the static 
gravitational field dn/dx = a>g. In the latter, u is the total energy of the test particle 
and g is the three-dimensional gravitational field strength vector. Although these equations 
are formally similar, there is a fundamental difference between them: the electromagnetic 
coupling constant e of a charged particle is not affected by its motion, while the coupling 
to the gravitational field is proportional to the total energy of the test particle. The latter 
property is a direct consequence of the equivalence principle. Let us formally replace the 
electric field strength E by a quantity that characterizes the gravitational field strength g 
and exploit the equivalence principle by considering, that the coupling of the particle to the 
field is proportional to the energy of the former, which also allows us to replace e by u. The 
expression that arises from ([33|) as a result of these replacements can be interpreted as the 
mean numbers of particles created by the corresponding gravitational field and have the form 
of the Boltzmann distribution with characteristic temperature 9' = If g is the gravita- 
tional field strength at the horizon of the black hole g(H), then 9' is only 2 time greater than 
the Hawking temperature 9(h)- I n spite of the fact that an explicit progress was achieved 
in the way of comparing both distributions, some questions remain. For example, why the 
temperature derived by means of the equivalence principle from the electrodynamical dis- 
tribution differs by the factor 2 from the Hawking one? Is there a thermal interpretation of 
the electrodynamical distribution (|33|) or some universal form of particles creation spectrum 
which is valid in both cases? Below we propose some interpretation of the electrodynamical 
formulas which pretends to answer these questions. We go beyond the classical considera- 
tion, taking into account properties of the physical vacuum in the time dependent external 
field. 

First, one can remark that due to the time dependence of the potential Ad(x°), which 
defines the quasiconstant electric field, the level of the vacuum energy changes with time. 
Thus, one has to calculate carefully the difference between the energies of the system in the 
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initial (vacuum) and in the final (with particles) states. Let us do the calculations in the 
case of fermions and in zero order with respect to radiative corrections in the T-constant 
electric field. In this case the corresponding Hamiltonian has the form 

H(x°) = J V(x)H . p .V(x)dx , (56) 

where H op was defined in (^), and ^f(x), \l/(x) are electron-positron field operators in the 
generalized Furry picture [p^9| , p0| ,|6l . Being written in terms of IN- and OUT- operators of 
creation and annihilation at a; — > =Foo respectively, the Hamiltonian H(x°) has the diagonal 
forms, 



H(x°) = ^Poih) a\ )r (in)a V)r (in) + b pr (in)b Pjr (in) - 1 
H(x°) = y^Pofa) a pr (out)a P;r (out) + b pr {out)b Pir (out) — 1 



x 



X ' 



-oo, 



+oo, 



(57) 



p,r 



where, as before, t 2 = —t\ = -|, and po{U) = y m 2 + p\ + (7c D (ti)) 2 is the energy of a particle 
in the initial and final time instants £j (the longitudinal momenta TToiti) = Pd — eAo(£i) in 
the T-constant field have the form ttd(±|-) = Pd T e ^>\- Let us consider the variation of 
the total energy of the system, which goes over from the initial vacuum state |0, in > to the 
final state \N, out > with pairs created in all the possible levels for the T-constant field, 

\N,out>= Y[ a pr (out)b pr (out)\0,out > . 

p±,r,\ Pn \<eET/2 

Then one can formally write the energy of the initial state as 



p.r 



and of the final state as 



2 Po (t 2 )-j2po(t2) 



\ PD \<eET/2 



PD 



Thus 



A£ = £ 2 - £ x = £ \Po(h) +Po(h)] + A£ v 

P^' r \\p D \<eET/2 



(58) 
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where 

A£ mc = E Ipoi^-poih)] (59) 

\ PD \>eET/2 

is the shift of the vacuum energy related to the levels with given p±,r, in which no pairs 
appear. We are going to analyze this shift only. That is why we do not discuss here 
regularization problems of the total sum (|58|) (that can be done, using, for example, the 



methods described in [111). One can see that 



E [Po(h) - Po(t 2 )] = o, 

\ PD \<eET/2 

since Po{ti) — Po(^) is an odd function of po- That allows one to extend the summation 
in fl5*P| ) over all the longitudinal momenta. The vacuum before the time instant ti was free 
and therefore symmetric with respect to the longitudinal kinetic momentum Kniti) = kd = 
Pd + eE^. Replacing the summation over pz> by one over 7Td, one can therefore treat the 
corresponding improper integral in sense of its principal value. Thus 

A£ mc = lim f M (^/ m 2 + P l + n 2 D - Jri* + pi + (vr D - eETf) d<K D 

= ~^{eETf = - A [n D (t 2 ) - n D (ti)} 2 ■ (60) 

Z7T Z7T 

Since the number of states with given p D , in which particles can be created, is equal to 
^eELT, see ([3^), then the shift ( |60|) can be rewritten in the form 

A£ mc = Ae„ ac , Ae mc = —eET = —\n D (t 2 ) — 7i£>(ti)|. (61) 

\ PD \<eET/2 

Thus, 

A£ = E E Ae ' Ae = Po(* 2 ) + Po(h) + Ae mc , (62) 

P^' r \PD\<eET/2 

where Ae can be interpreted as a work which the external field accomplishes for the creation 
of a pair in a given state. It contains a contribution Ae mc which takes into account a shift 
of the vacuum energy in those states which remain vacuum ones. The corresponding work 
with respect to a particle will be denoted by u, so that, 
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1 . 

1 

~ 2 



m 2 + pi + (n D (t 2 )) 2 + Jm 2 + p 2 ± + (Mh)) 2 - - Mh) 



(63) 



Now we remark that due to the conditions of the stabilization T >> To, \pd\ < e -^f"> under 
which all the results for the T-constant field were obtained, we can write 

1 / 1 1 \ A XeE m 2 + p 2 ± 

u = -XeE - — — + — — = - = --, A = . 64 

4 \\M*2)\ \Mti)\) T 2p (t 2 y eE 1 ; 

Then the spectrum (33) can be rewritten in the following form [] 



UJ 



Np, r = exp |— 27r-|— | , (65) 

where the quantity g can be written in several equivalent forms 

_ceE ( 1 1 \ _ 2c _ ceE 



The last expression in ( pEp allows one to treat g as the classical acceleration of a particle 
in the electric field in the final time instant t 2 = for the case when the action time 
of the field is big enough, so that the corresponding velocities are near c. Formally this 
is valid under the quantum condition of stabilization. The distribution ( |65|) is, in fact, 
the Boltzmann one with the temperature 9 = 2 nck B navm § literally the Hawking form. 
Thus, if one identifies the work u, we have introduced, with the energy of a particle in the 
formula (|55|), then the distributions in electrodynamical and gravitational cases have the 
same thermal structure. Let us discuss now the possible origin of the differences in the 
electrodynamical and gravitational formulas. First of all, the formula fl55D is derived in the 
formalism of the stationary scattering theory, where it is not necessary to take separately 
into account a shift of the vacuum level. In this case the energy of a particle created may 
coincide with the corresponding work of the field. Second, the different form of the thermal 
distributions (Boltzmann, Planck) can be stipulated by essentially different situations in 



L we have restored h and c here for convenience of the reader 
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both cases. In the electrodynamical case we deal in fact with pure states, whereas in the 
gravitational problems a density matrix is arisen necessary due to the horizon of events 
formation. At u/g << 1 the Planck spectrum coincides with the Boltzmann one. In this 
case one can believe the form ( |65| ) for the spectrum of particles created is universal and 
applicable to any theory with quasiconstant external fields. 

The form (^) can be useful to describe situations in constant fields, where it is convenient 
to avoid the consideration of the time evolution, as we have seen comparing it with the 
gravitational cases. Another form of the distribution (^) with 2u = Ae from (^) (or 
(|62|)) and with the acceleration g from (|66D, can be useful in problems with explicit time 
dependence. 

The universality of the formula (^) can be examine also in the case of the adiabatic 
electric field, a » a , considered in Sect.IIIB. To apply it to the latter case one needs to 
put ^2,1 — * ±00, then 7T£i(t 2 ) = pn — eEa, 7Td(£i) = Pd + eEa and pofe) = w+, Po(ti) = 
In this case according to ( |61~D Ae mc = — 2eEa and the distribution ( [46]) follows. 

One can remark in this connection that in case E/E c << 1 («o << T ) the formulas (|B3|) 
(and, therefore, (§^)) in the adiabatic field are valid at the condition t 2 = —ti = ^ >> «o ~ 
Tjf , which are weaker than in the case of the T-constant field. That allows one to interpret 
the time Tjf , which already had appeared before in Subsect.A as a pair formation time, 
from another point of view. Considering the formula (|63D , one can see that at T >> T/ 
the Ae is less than 2^Jm 2 + p 2 , due to the essential contribution of the vacuum shift Ae vac . 
That means the work of the external field to produce a pair is less than one, which could be 
expected from the perturbation theory, where no vacuum change is taken into account. 

The consideration presented was made only for fermions. However, if one believes that 
the quantity Ae mc can be taken in the form ( |6T]) for bosons as well, then the distribution ( |65| ) 
holds also in the scalar case. A consistent analysis for charged boson is more complicated 
and needs to take into account possible condensate formation and its evolution in an external 
field (see, for example, PJ4"0|]). 

Finally, we believe that the formulas derived and the time-dimensional analysis presented 
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can be also useful to describe some collective effects in the framework of quantum field theory, 
for instance, to describe multiple particles creation by means of the external field approach 
T3| , [T4|| or in string models with external field |TTyi2]j4~If . 
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